Anomalous Transport and Possible Phase Transition in Palladium Nanojunctions 
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Many phenomena in condensed matter are thouglit to result from competition between different 
ordered phases. Palladium is a paramagnetic metal close to both ferromagnetism and supercon- 
ductivity, and is therefore a potentially interesting material to consider. Nanoscale structuring 
of matter can modify relevant physical energy scales leading to effects such as locally modified 
magnetic interactions. We present transport measurements in electromigrated palladium break 
junction devices showing the emergence at low temperatures of anomalous sharp features in 
the differential conductance. These features appear symmetrically in applied bias and exhibit 
a temperature dependence of their characteristic voltages reminiscent of a mean field phase 
transition. The systematic variation of these voltages with zero-bias conductance, together with 
density functional theory calculations illustrating the relationship between the magnetization of 
Pd and atomic coordination, suggest that the features may result from the onset of spontaneous 
magnetization in the nanojunction electrodes. We propose that the characteristic conductance 
features are related to inelastic tunneling involving magnetic excitations. 



Introduction 

Many phenomena in modern condensed matter 
physics, including metal-insulator transitions, high tem- 
perature superconductivity, heavy fermions, multiferroic 
effects, and quantum criticality, are thought to orig- 
inate from competition between different possible or- 
dered electronic ground states. Slightly shifting the bal- 
ance between these orders via external fields, strain, dis- 
order, or carrier concentration can tune a system be- 
tween different phases. Palladium's proximity to both 
ferromagnetism^"— and superconductivit}*^"— make it an 
interesting material to consider in this context. Structur- 
ing matter on the nanometer scale is one means of modi- 
fying relevant physical energy scales, with nanoscale con- 
finement already known to favor locally modified mag- 
netic interactions in Pdi^i^rJ^ Nanogap structures can be 
particularly useful in examining the relevant physics in 
such complicated systems. Tunneling transport has been 
a valuable spectroscopic tool, with the tunneling conduc- 
tance as a function of bias voltage typically interpreted 
as probing the convolution of the electronic densities of 
states (DOS) of the two electrodes. In fixed nanogap 
structures, this has given insight into varied systems, in- 
cluding charge density wave compoundssii and assorted 
superconductors J^iH 

In this article we report tunneling measurements in 
Pd-Pd tunnel junctions fabricated by electromigration 
of planar constrictions on oxidized Si substrates. We 
find a characteristic pattern of features in the differential 
conductance appearing at source-drain voltages {Vsd) lin- 
early related to the log of the zero-bias conductance (ZBC 
= G{Vsd = 0)). These features rapidly decrease in bias 
with increasing temperature, with a temperature depen- 
dence reminiscent of a mean-field transition. Of particu- 
lar interest is the contrast in scale between the energy of 



an individual peak or dip (eV^d) and the respective tem- 
perature at which it vanishes (fc^Tc), which often differ 
by ^ 100 X. These features are robust and appear to be 
particular to Pd electrodes. We discuss possible physical 
mechanisms for these features, and suggest that they are 
likely signatures of inelastic tunneling processes associ- 
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FIG. 1: (a) Scanning electron microscope (SEM) image of 
a Pd nanowire on Si02. Inset: High resolution image af- 
ter electromigration exhibiting a tunneling gap separation < 
2 nm and an absence of nanoparticles. (b) dl /dV vs. Vsd, 
at T — 1.8 K, for representative sample 1, which includes a 
pair of broad high bias peaks near ±215 mV and a series of 
narrow features at lower bias, nearly symmetric about l/gd = 
V. (c) Waterfall plot of the lower bias region of sample 1, 
showing the temperature evolution. Temperature is increased 
from 1.8 K (bottom) to 4 K (top) in steps of 100 mK. Traces 
are offset for clarity by 0.006 e /h. 



ated with the onset of spontaneous magnetization in the 
electrode tips. While the details have not been definitely 
determined, calculations of the expected magnetic mo- 
ment in model Pd nanocontacts are consistent with both 
the onset temperature of the conductance features and 
the observed systematic variation of voltage scales with 
the ZBC. 

Results 

Bowtic-shaped Pd nanowires are fabricated on a 
Si/Si02 substrate (see Methods). The devices are 
cleaned and subsequently cooled to < 2 K. A tunneling 
gap is created in a nanowire by employing the electromi- 
gration technique (Fig. la).— Differential conductance 
(dl/dV) is measured as a function of V^d, and the con- 
ducting substrate is used as a back gate by applying a 
separate bias, Vg. 

When the electromigration process reduces the ZBC 
below Go (most commonly between O.SOGo and O.OOIGq, 
where Go = e'^/h), a characteristic series of peaks and 
dips is observed in the curve of dl/dV vs. V^d when 
measured at low temperatures. These differential con- 
ductance features (DCFs) are highly reproducible in Pd. 
Out of 111 clean Pd nanowire samples successfully elec- 
tromigrated to a measurable ZBC less than Gq, more 
than 75% of devices exhibited characteristic DCFs and 
their corresponding temperature dependence. This holds 
true for samples fabricated on both n- and p-type Si sub- 
strates. When the electromigrated nanojunctions have 
a ZBC of G{Vsd = 0) > Gq, these features are entirely 
absent. These DCFs appear to be specific to Pd. Ex- 
periments with both Ni-Ni and Au-Au tunnel junctions 
formed via analogous fabrication processes have shown 
no indications of similar DCFs. 

When T ~ 1.8 K we distinguish distinct types of DCFs 
associated with two different bias regimes. In the "low 
bias" region (typically |V^d| ^ HO mV), a series of sharp 
features appears on each side of V^d = V (Fig. lb & c) . 
From device to device, these peaks and dips occur with 
varying amplitudes, lineshapes, and spacings; for a given 
device, the features are clean and reproducible. The 
"high bias" region (typically jV^dl ~ 150 mV - 400 mV) 
contains the second type of feature, which is always man- 
ifested as a single pair of relatively broad peaks, with one 
peak appearing on each side of 14d = V (Fig. lb). The 
amplitude of these two peaks is greater than those in 
the low bias region. While the pattern of all DCFs ap- 
pear symmetric about T4d = V, there can be an overall 
asymmetry in conductance magnitudes due to a back- 
ground contribution (roughly linear in V^d) to dl/dV. 
Such a background is commonly observed in STM and 
break junction experiments J^i^ 

The DCFs undergo a rapid decrease in voltage position 
and amplitude as temperature is increased (Fig. Ic). At 
sufficiently high temperatures all temperature-evolving 
features have disappeared, and the resulting trace of 
dl/dV appears weakly nonlinear, consistent with what is 
expected for an ordinary metallic break junction device 
with a small interelectrode gap. A critical temperature. 
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FIG. 2: Base temperature traces of dl /dV vs. VLd and their 
corresponding temperature dependence, (a) dl/dV vs. Kd 
for device 2. (b) Colormap depicting temperature evolution 
of device 2. Conductance amplitude is indicated according to 
colorbar on the right. The high bias peaks terminate near 9 K. 
The lower bias features evolve in a similar fashion, but are not 
observable on this color scale, (c) dl /dV vs. T4d for device 
3 showing only the features in the low bias regime. Many 
peaks and dips are evident, (d) Colormap of temperature 
evolution for device 3. Each feature in the spectrum of dl /dV 
disappears at a different value of T, which we refer to as its 
critical temperature, Tc. 



Tc, can be assigned to each DCF, defined as the temper- 
ature at which that particular feature approaches VJsd = 
V. The energy of an individual peak or dip cannot usu- 
ally be tracked all the way to eV^d = V as it invariably 
becomes difficult to resolve with diminishing amplitude, 
resulting in the apparent "smearing" of the colormaps in 
the vicinity of the Kid = V axis (Figs. 2, 3, and 4). 
The energies of the DCFs decrease with temperature in 
a manner that is well fit by the equation: 



eVsdiT) = A 1 - — 



T 
Tc 



(1) 



where A is the voltage position of the peaks (multiplied 
by e) as T ^- K and /3 = i, the classical exponent for 
a mean-field transition in a Landau theory. The low bias 
DCFs of all Pd devices have T^ < 4.3 K and evolve with 
temperature as in Equation ([T|). An example is shown 
in Fig. 3, along with fits to Equation ([T|). The energy of 
the high bias peaks also follows this form, but begins to 
deviate at T/Tc < 5, where T^ is typically between 8-12 
K (Fig. 2b and Supplementary Fig. S5). 

In Figure 4 the voltages of representative peaks and 
dips, from a set of several dozen devices, are plotted as 
a function of the log of their corresponding ZBC when 
T « 1.65 K. The peaks associated with high bias regions 
are quite linear in \og{G{Vsd = 0)), indicating that the 
voltage position of these features decays exponentially 
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FIG. 3: Temperature dependence of the voltages at which DCFs appear is well characterized by equation ([T]). (a) dl /dV vs. 
\4d at T = 1.65 K for sample 4. (b) Colormap of sample 4 formed from traces of dl /dV vs. Kd acquired at temperatures 
between 1.65 K and 3.325 K in steps of 25 mK. (c) Plot of voltage at which the two primary conductance peaks (on both sides 
of Vsd = V) appear versus temperature. Solid lines are least squares fits to the data using equation ([T]). Error bars represent 
uncertainty in peak position of Gaussian lineshape fit to each peak using a least squares method. 



with ZBC. Fits to the high bias data indicate that these 
peaks will go to zero bias as G{Vsd = 0) — )- Gq. This 
is consistent with our observation that DCFs are only 
present for samples with G{Vsd = 0) < Gq. DCFs in the 
plot associated with the low bias regime appear anywhere 
within a bias range between Vsd = V and an upper |T4d| 
Umit. The edge of this envelope is also roughly linear in 
the log of G{Vsd = 0). 

The ZBC of a tunnel junction with a vacuum gap is 
expected to decrease exponentially with increasing inter- 
electrode gap distance, d, according to G a e~°"^, where 
a oc y/(j), and ^ is a typical barrier height related to the 
metal work function. Since 14d oc log(G'(ysd = OT^)), a 
linear relationship exists between the 14d positions of the 
DCFs and d. This is further confirmed by observations in 
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FIG. 4: Exponential relationship between voltage position of 
DCFs and corresponding ZBC. Vsd position of representative 
dl /dV features at T ~ 1.65 K for many junctions plotted as 
a function of the log of G{Vsd = 0) for each sample. Points 
pertaining to high bias peaks follow linear curve. Black line 
is a linear fit to positions of the high bias peaks for Vsd > 
V. Points pertaining to low bias features appear between I4d 
= V and an envelope outlined by an approximately linear 
curve. 



a single device when measurements of dl/dV vs. 14d are 
performed after successive electromigrations: The ZBC 
is reduced and the DCFs are found at higher biases while 
the pattern of DCFs remains relatively unchanged (see 
Fig. S2 in Supplementary Information). 

Conversion of the ZBC data in Fig. 4 into a dis- 
tance scale is extremely sensitive to both the conduc- 
tance prefactor and a, neither of which is well estab- 
lished in Pd.— If desired, one could perform theoretical 
calculations based on density functional theory to obtain 
a quantitative estimate for a with particular tip geome- 
tries. A typical value of a for vacuum tunneling between 
high work function metals is 2 A~^^ Using this value 
and assuming a prefactor for the conductance of 2e^ /h 
when the two electrodes are in contact at the single atom 
level, we then estimate that d ranges between 1 and 4 A. 

We have also examined the effects of gate bias and 
magnetic field on the junctions, as further explained in 
the Supplementary Information. As Vg is increased in 
the positive (negative) direction, the ZBC decreases (in- 
creases), and the DCFs move to higher (lower) energies, 
increasing (decreasing) their associated Tc's (Fig. Sic). 
Some DCFs appear to split at large negative Vg. There 
is no sign of Coulomb stability diamonds, and these gate 
bias effects have the same polarity for both n- and p-type 
substrates. External magnetic fields up to ±9 T in a di- 
rection perpendicular to the plane of the devices lead to 
no significant effects. There is no Zeeman splitting or 
broadening of peaks in dl /dV , and no hysteresis seen in 
measurements of magnetoresistance (Fig. Sid). There 
is a very slight shift (~ ImV) of DCFs to higher volt- 
ages at ±9 T, but we cannot rule out that this is due 
to slight field-dependent changes in sample temperature. 
SQUID measurements of co-deposited Pd films show only 
paramagnetic response down to 2.0 K. 

Discussion 

Selectively altering matter on the nanometer scale is 
one means of modifying relevant physical energy scales 
and can lead to collective electronic effects not necessar- 



ily present in bulk samples. The anomalous transport 
data presented here, found in a chemically homogeneous 
system, are an example of such an effect. We consider 
several possible explanations for the observed phenom- 
ena. One concern is that the DCFs are extrinsic rather 
than an inherent property of the Pd junctions. For exam- 
ple, the presence of helium exchange gas in the chamber 
raises the possibility that an adsorbed He film could be 
the origin of the features. This is unlikely however, since 
the DCFs are also observed when no He exchange gas is 
used. Also, with or without exchange gas, the high bias 
DCFs survive well above 4.2 K. 

Similarly, it is well established that hydrogen has an 
affinity for PdJ^ Previous experiments involving tunnel- 
ing spectroscopy of metallic nanogaps in the presence of 
H2 reveals low temperature dl/dV features reminiscent 
of our data due to transitions between coverage depen- 
dent statesi^ or vibrational modea^S of the molecules. 
The voltages of the DCFs we report arc similar in scale 
to the energy of phonons or vibrational modes related to 
adsorbed molecules. However, phonon mode energies do 
not have the strong temperature dependence exhibited 
by these DCFs. Furthermore, unlike the H2-containing 
junctions, our devices are stable over time (3-4 days), 
do not exhibit enhanced noise or telegraph switching, 
and the DCFs do not appear with highly reproducible 
energies 1^ 

One possible interpretation of the data is that dl/dV 
reflects a convolution of the (thermally broadened) DOS 
of the source and drain electrodes or a particle between 
the electrodes, as in scanning tunneling spectroscopy. In 
this situation the DCFs would correspond to features in 
the DOS at energies eKid away from the Fermi level (Ep). 
A probe of the electrode DOS would not be expected 
to produce symmetric features with respect to Vgd, nor 
would there be reason for the peaks to shift towards Ep 
with increasing T. The symmetry of dl/dV features with 
respect to 14d = 0, together with the observed gate and 
temperature dependencies are also inconsistent with ex- 
pectations of Coulomb blockade physics {e.g., resonant 
tunneling through a nanoparticle in the interelectrode 
gap). 

Pd can undergo a superconducting phase transition 
when doped with H^ or when irradiated with He ionsi 
to create disorder and lattice strain. The form of the T- 
dependence of the DCFs in our devices is the same as that 
used to describe the energy gap, A, in superconducting 
materials. The dl/dV data have some qualitative resem- 
blance to that seen in superconductor-superconductor 
tunnel junctions. However, as temperature is increased 
one would expect all dl/dV peaks to approach a common 
value of Tc, contrary to the behavior of the DCFs in our 
Pd break junctions. Moreover, the gap energies inferred 
from interpreting the DCFs as DOS features are orders 
of magnitude away from what one would anticipate given 
the observed TcS. 

Figure 4 imposes a significant constraint on any candi- 
date explanation for the data. It is difficult to construct 
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FIG. 5: DFT calculations of magnetic moments for the end 
three layers forming the tip of a Pd electrode grown in the 001 
direction as a function of the tip-tip distance, d. Figures in 
the inset represent the model tips for contact (upper, d = A) 
and tunneling (lower) regimes. Results show that Pd tips are 
magnetic, but with atomic magnetic moments smaller than 
those found in Pd chains (for the same type of calculation) 
due to higher coordination of the atoms. 



a scenario in which the energies of features in the DOS 
vary systematically with interelectrode conductance as 
observed. In a typical tunnel junction, increasing the 
electrode separation will reduce the conductance ampli- 
tude related to features in the DOS {e.g. from resonant 
tunneling into accessible states), but should not alter 
their respective energies. Fig. 4 strongly suggests that 
the DCFs are not simple DOS resonances. 

The linear relationship between DCF Vsd and the in- 
ferred interelectrode separation d suggests the possibility 
that the DCFs represent physics driven by the interelec- 
trode electric field, Vsd/d. However, any electric field- 
based mechanism is constrained by the mean field-like 
temperature dependence seen in Fig. 3. 

One candidate mechanism would be electric field con- 
trol of magnetic ordering, which has been suggested by 
several groups22i2ii^ and demonstrated experimentally 
in FcPdi^ Indeed, hole-doping of Pd with an electric 
field has been predicted to instigate a transition to fer- 
romagnetism, causing a step-like change in the tunneling 
current at a given value of applied bias.— In our struc- 
tures, the bias across the tunnel junction induces a local 
(on the scale of the screening length) accumulation of 
electrons and holes on the facing electrodes, resulting in 
a raising or lowering of the Fermi level at the surface, re- 
spectively. Palladium has a peak in the DOS just below 
the Fermi level. It is conceivable that a sufficient accu- 
mulation of holes could result in the satisfaction of the 
Stoner Criterion leading to itinerant fcrromagnctism at 
that surface. 

However, there are a number of problems with this 
scenario. For example, this idea fails to account for the 
symmetry observed in the conductance spectrum about 
14d = 01^. Which electrode is hole-accumulating depends 



on the bias polarity. It seems unlikely that the two po- 
larities would result in identical electric field-driven tran- 
sitions. It is also unclear how such a mechanism would 
lead to the rich variety of DCFs and their observed de- 
pendence on Vg (See Supplementary Information). Fur- 
thermore, calculations indicate that the strength of the 
field necessary to instigate such a transition (> 1 V/A) 
would require an applied bias prohibitively high in our 
devices i^i^^ Therefore, we do not think that electric field- 
driven magnetism is a likely explanation for our observa- 
tions. 

We speculate that the observed DCFs are related to 
inelastic electron tunneling. An inelastic tunneling ori- 
gin for the low bias DCFs would naturally explain the 
observed symmetry in Vgd , and the variety of lineshapes 
could be explained by interference between various elas- 
tic and inelastic channels involving local modes.— i^^ In 
contrast, the high bias DCFs are exclusively broad max- 
ima, suggesting that different processes may be involved 
for low and high bias features. However, the fact that 
both regimes are strongly correlated with one another 
(Fig. 4) and have the same trend with zero-bias conduc- 
tion strongly implies that the low and high bias DCFs 
share crucial underlying physics. 

We further propose that the relevant transport chan- 
nels may arise due to the onset of spontaneous magne- 
tization at the electrode tips due to the proximity of 
Pd to a magnetic instability. While Pd surfaces are not 
magnetic for film thicknesses exceeding a few monolay- 
ers, reports suggest the possibility of itinerant ferromag- 
netism in Pd nanoparticles^''^ and atomic-size contacts^ 
Zero-tcmpcrature density functional theory (DFT) calcu- 
lations for model Pd nanocontacts, as those shown in Fig. 
5, indicate the possible presence of magnetism due to the 
low coordination of the atoms forming the electrode tips. 
(See Supplementary Information for calculation details). 
A quantitative estimate of the Curie temperature is, how- 
ever, beyond DFT capabilities. Furthermore, even if we 
take the calculated energy difference per atom between 
the ferromagnetic and paramagnetic states as a reference 
value, the inherent uncertainties in the calculations are 
larger than the range of the observed DCF TcS. Nev- 
ertheless, we can say that this energy difference sets an 
upper limit to Tc that is consistent with the observed 
values (sec Supplementary Information). 

If one assumes an excitation energy proportional to the 
local magnetization (consistent with, e.g., Stoner excita- 
tions), then the temperature evolution of eV^d would cor- 
respond to the onset of magnetic order. Moreover, within 
that assumption the same model calculations are consis- 
tent with the trend of Fig. 4. The layer-by-laycr mag- 
netic moment per atom (averaged over layer), < (si) >, 
is plotted in Fig. 5 for the three atomic end layers of 
one tip (the other tip is identical) as a function of the 
distance between tips, d. The calculated average mag- 
netic moment per atom decreases as one moves away from 
the ends of the tips towards the bulk, revealing the ex- 
pected inhomogeneous nature of the magnetization. The 



magnetic moment of the tip atom is visibly affected by 
the interelectrode distance, decreasing with d below 5 A, 
consistent with the dynamic range of d estimated ear- 
lier from the measured ZBC. This can be rationalized by 
noting that on increasing the coordination (decreasing d) 
the DOS decreases at the Fermi level and the polariza- 
tion becomes weaker, in accord with the Stoner criterion. 
The electrodes in a real break junction may contain many 
undercoordinated atoms, giving rise to localized areas of 
increased magnetic susceptibility, extending over a larger 
region of the junction. In this picture, DCFs with differ- 
ent TcS would result from different local portions of the 
electrodes. 

Another related speculative mechanism for DCFs in- 
volves the interplay of spontaneous magnetization, as 
above, and Kondo physics. The formation of the highly 
correlated Kondo state in bare metal break junction 
devices22i^ can lead to sharp features in dl/dV at sym- 
metric voltages. In this case the dl/dV peaks would cor- 
respond to a spin-split zero-bias resonance resulting from 
a strong local exchange field established by the presence 
of ferromagnetisniiSi^ or effective magnetic impuritiesi^ 
The T-dependent onset of local magnetization at the elec- 
trode tips could lead to such a splitting that increases 
with decreasing temperature. 

The apparent lack of magnetic hysteresis and signifi- 
cant magnetoresistive effects would seem to cast doubt on 
a hypothesis involving magnetic correlations. However, it 
is worth recalling that Stoner-like excitations would not 
Zeeman split, although they may be expected to shift 
with the application of an external magnetic field. Fur- 
thermore, magnetically ordered regions possibly present 
very small local anisotropy energies^ and a superparam- 
agnetic response to externally applied fields. This implies 
a vanishingly small coercivity and suggests that large hys- 
teretic effects are not necessarily expected.— The mis- 
match between DCF energy scale and Tc is not readily 
explained, and would require a more detailed theoretical 
picture of the inferred magnetic transition and result- 
ing excitations. One should recall, however, that Stoner 
excitations (invoked as a possible explanation of DCFs) 
and the spin stiffness (related to Tc) are associated with 
different exchange couplings ii2^ This makes their charac- 
teristic energy scales differ by an order of magnitude, at 
least in transition metals. We see no reason to expect 
otherwise in the case of Pd. 

The data presented here unambiguously show that 
unusual electronic effects are present in Pd nanojunc- 
tions at low temperatures. The trends in the data with 
T, G{Vsd = 0), Vsd, and Vg greatly restrict possible 
mechanisms. While speculative, an explanation of the 
features in terms of inelastic processes involving mag- 
netism is consistent with the observed temperature range, 
source-drain symmetry, special character of Pd, and the 
trend with inferred electrode separation. This idea is 
also testable through doping the Pd with a ferromag- 
netic metal to see if enhancing the tendency toward fer- 
romagnetism alters the phenomenon. These conductance 



anomalies demonstrate once again that even "simple" 
chemically homogeneous materials can exhibit rich and 
surprising phenomena at the nanoscale. 

Methods Arrays of bowtie-shaped nanowirc patterns 
are defined by electron beam lithography with minimum 
widths between 80 and 100 nni. Samples are composed 
of 14 nm of Pd deposited by c-bcam evaporation on a 
highly doped Si substrate with a 200 nm Si02 insulating 
layer (Fig. la). The devices are contacted individually 
using a piezo-controUed Attocube^'^'^probe station inside 
of a ''He cryostat. An array of devices is cleaned in sol- 
vents, exposed to an oxygen or argon plasma clean, and 
subsequently inserted into the sample space of the probe 
station, which contains approximately 10 mB of He ex- 
change gas at 300 K. The electromigration technique^^ 
is employed at low temperature (T < 2 K) to create a 
small gap in the nanowire (Inset, Fig. la). Differential 
conductance {dl /dV) is measured as a function of source- 
drain bias (Kid) using standard low frequency lock-in 
techniques with an rms excitation voltage of 625 /iV. T 



V is also measured using a voltage source and current 
preamplifier, for a small portion of samples. 
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